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Abstract. We discuss the single transverse-spin asymmetry (SSA) to be observed in the D-meson 

production with large transverse-momentum in semi-inclusive DIS (SIDIS), ep^ — s> eDX. This 

contribution is embodied as a twist-3 mechanism in the collinear factorization, which is induced 

by purely gluonic correlation inside the transversely-polarized nucleon, in particular, by the three- 

gluon correlation effects. We derive the first complete formula for the corresponding SSA in the 

"^ ' leading-order QCD, revealing the five independent structures with respect to the dependence on the 

azimuthal angle for the produced D-meson. Our result obeys universal structure behind the SSAs 

O ■ in a variety of hard processes. We present the numerical calculations of our SSA formula for the 

QJ I D-meson production at the kinematics relevant to a future Electron Ion Collider. 
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>, 
'nI" ! Charm productions in semi-inclusive DIS (SIDIS), as well as in the pp collision, are 

known to be associated with the twist-2 gluon distributions in the nucleon, since the cc- 
cr] • pair creation through the photon-gluon or gluon-gluon fusion is their driving subprocess. 

(^ . Similarly, the twist-3 contributions in the charm productions can be generated by the 

O I purely gluonic effects inside the nucleon, in particular, the multi-gluon correlations. 

Indeed, the observation of the single transverse-spin asymmetry (SSA) in the open charm 
productions allows us to probe the corresponding twist-3 effects [1, 2, 3, 4, 5]. 

The corresponding SSA arises as a naively T-odd effect in the cross section for the 
/\ • scattering of transversely-polarized nucleon off an unpolarized particle, observing a D- 

j^ ■ meson with momentum Ph in the final state, and this requires, (i) nonzero transverse- 

momentum P/,^ originating from transverse motion of quark or gluon; (ii) nucleon he- 
licity flip in the cut diagrams for the cross section, corresponding to the transverse polar- 
ization; and (iii) interaction beyond Born level to produce the interfering phase between 
the LHS and the RHS of the cut in those diagrams. In particular, for large Pj^^ ^ Aqcd 
to be dealt with the collinear- factorization framework, (i) should come from perturbative 
mechanism as the recoil from the hard (unobserved) final-state partons, while nonpertur- 
bative effects can participate in the other two, (ii) and (iii), allowing us to obtain observ- 
ably large SSA. This is realized with the twist-3 multi-gluon correlation functions for 
the transversely-polarized nucleon [3]. This twist-3 mechanism may be considered as an 
extension of the corresponding mechanism for the SSA in the pion productions in the 
SIDIS [6], pp collisions [7], etc., based on the quark-gluon correlations in the nucleon, 
but it has been clarified [3] that a straightforward extension [1,2] leads to missing many 



FIGURE 1. Feynman diagrams for partonic subprocess in ep^ — > eDX\ mirror diagrams also contribute. 



terms in the SSA. We shall discuss the leading-order (LO) QCD result for the SSA in the 
high-P/.j, D-meson production in SIDIS, e(£) +p{p, S^) -)■ e(f ) +D{Pi,) +X, controlled 
by the twist-3 three-gluon correlation functions in the nucleon. We also present a numer- 
ical estimate of the corresponding SSA with kinematics of a future Electron Ion Collider 
(EIC). We use, as usual, the kinematic variables Sep = {i + p)^, q = i — i', Q^ = —q^, 
Xhj — Q^/{2p ■ q), and Zf = p- Ph/{p- q)- We work in a frame where the momenta q and 
p are coUinear, and define qr = Ph±/zf and the azimuthal angles ^, $5, and x of the 
lepton plane, the spin vector S^, and the D-meson momentum PJ^, respectively [3, 4]. 
We take into account the masses m^ and m/^ for the charm quark and the D meson. 

In the LO in QCD perturbation theory, the photon-gluon fusion subprocesses of Fig. 1 
drive the SSA, where the above-mentioned (i) is provided by the recoil from the hard 
unobserved c quark and the c quark with the momentum pc fragments into the D-meson 
in the final state. The short bar on the internal c-quark line indicates that the pole part 
is to be taken from that propagator, to produce the interfering phase of (iii); we note 
that these pole contributions from Fig. 1 would cancel the similar contributions from 
the corresponding mirror diagrams, if the c quark were unobserved in the final state 
as in the D-meson production. The external curly lines represent the gluons that are 
generated from the three-gluon correlations present inside the transversely-polarized 
nucleon, {pS±\Aa{0)Ap{ri)AY{^)\pS±), corresponding to (ii). The diagrams obtained 
by the permutation of the gluon lines in Fig. 1 also satisfy (i)-(iii), but the Bose statistics 
of the gluons in the above matrix element guarantees that we need not consider those 
diagrams separately. Thus, the SSA in the present context can be derived entirely as the 
contributions of soft-gluon-pole (SGP) type [6], leading to k2 — ki = 0, by evaluating 
the pole part in Fig. 1 . The twist-3 nature of those contributions are unraveled by the 
coUinear expansion, as usual. The expansion produces lots of terms, each of which is not 
gauge invariant. Indeed, many of them vanish or cancel eventually, and the remaining 
terms can be organized into a gauge-invariant form. This can be demonstrated [3] by 
sophisticated use of the Ward identities for the contributions of the diagrams in Fig. 1. 
The resulting factorization formula of the spin-dependent, differential cross section for 
ep^ —^ eDX reads [3, 4] ([dco] = dxbjdQ^dzfdqjd^dx, for the differential elements) 
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where x = x^j/x, z = Zf/z, and a complete set of twist-3 gluonic correlation functions 
is defined through the gauge-invariant lightcone correlation of three field-strength ten- 
sors [3] (rp- = 0, p ■ n = 1), 
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and similarly for A'^(jci,jc2) with the replacement J**^'" — )■ i/^*^^, and Dc{z) denotes the 
usual twist-2 fragmentation function for a c-quark to become the D-meson. The quark- 
flavor index a can, in principle, be c and c, with 5c = 1 and 5c — —1, so that the 
cross section for the D-meson production ep^ — )■ eDX can be obtained by a simple 
replacement of the fragmentation function to that for the D meson, Da{z) — )■ Da{z). The 
summation J^^ implies that the subscript k runs over 1 , 2, 3,4, 8,9, with =g/i = 1 H-cosh^ i//, 
^2 = —2, £^3 = — cos^sinh2i//, ^ = cos2^sinh i//, £^ = — sin^sinh2i//, and £/g = 
sin 2^ sinh^ \]/, where coshi// = IxyjSep/Q^ — 1 and ^ = ^ — X^ ^iid y^ is defined as 
yk = sin($5 - x) for fc = 1, 2, 3,4 and y^ = cos($5 - x) for k = %,9. Uem is the fine- 
structure constant, and Cc = 2/3 is the electric charge of the c-quark. Partonic hard parts 
Ad^ depend on nic as well as other partonic variables; for the explicit formulae of Ad^, 
we refer the readers to [3]. Note that, instead of evaluating the diagrams in Fig. 1 as 
above, those results can be obtained using the "master formula" [4], schematically given 
by 

d^Ao . fdx fdzdje^y{xp,q,pc) Aiva^ N^n^^ ^-2^ 

— ®Jip {x,x)(g)Da{z), (3) 
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where Jiff^v{xp,q,Pc) denotes the partonic hard part for the 2 — )■ 2 Bom subprocess, 
expressed by the diagrams in Fig. 1 with the soft (^2 — ^1=0) gluon line removed. This 
reveals that Ad^ in (1) are related to the twist-2 hard parts J^iv{xp,q,pc), similarly as in 
the SSA in various processes associated with twist-3 quark-gluon correlation functions. 
Reexpressing as ^ — x = ^ = ^h, ^s — X = ^h — ^s ^^ the above-mentioned explicit 
forms of ^^4 and ^^, where ^^ and 05 represent the azimuthal angles of the hadron plane 
and the nucleon's spin vector S±, respectively, measured from the lepton plane, one may 
express (1) as the superposition of five sine modulations. 
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with the corresponding structure functions /i , /2, • • • , /s, similarly as the twist-3 SSA for 
ep^ — )> enX, generated from the quark-gluon correlation functions [8]. 

As an illustration, we evaluate the SSA for the D^ production, ep^ — )> eO'^X, using 
(1), (4). In particular, we present Fi = /i/cf^ with the kinematics relevant to a future 
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FIGURE 2. The SSA as a function of (a) Zf and (b) Pt = Ph± at EIC kinematics with Sep = 2500 GeV^ 



EIC in Fig. 2, where o\^ denotes the twist-2 unpolarized cross section for ep — > eD^X, 
J^C7™P°'/[(ia)], averaged and integrated over the azimuthal angles 0/^ and x^ respec- 
tively. The solid and dashed curves show the contributions generated by the three-gluon 
correlation functions of (2), assuming 0{x,x) = O{x,0) = 0.005xG{x) [5] (see also 
[1, 2]) with the unpolarized gluon-density distribution G{x) for the nucleon and us- 
ing CTEQ6L parton distributions and KKKS fragmentation functions [9] with the scale 
jU^ = Q^ + rn^ +ziq-^. The results demonstrate good chance to access multi-gluon effects 
at EIC. The detailed numerical studies, including those for the asymmetries generated 
by the other structures /2, • ■ • ,/5 in (4), will be presented elsewhere. 
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